Theory of the NMR 1/Ti relaxation rate in a quantum spin nematic 
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The idea that a quantum system could behave hke a hquid crystal holds a longstanding fascination. 
One widely-discussed example is the quantum spin-nematic, in which magnetic moments order 
as spin-quadrupoles, breaking spin-rotation symmetry without breaking time-reversal symmetry. 
However, despite a long history, experimental evidence for the existence of the spin-nematic state 
remains scarce. Here we show how the NMR l/Ti relaxation rate can be used to distinguish the 
spin-nematic state from more conventional forms of magnetic order. Building on a symmetry-based 
theory of spin excitations, we find that, in a spin-nematic phase, l/Ti is dramatically suppressed, 
does not diverge at continuous phase transitions, and exhibits unconventional power laws at low 
temperatures. We discuss how this theory could be applied to LiCuV04, where spin-nematic order 
has been proposed to exist at high magnetic field. 

PACS numbers: 75.10.Jm, 75.40.Gb 
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The search for new, quantum phases of matter is cen- 
tral to modern condensed matter physics. Each new 
phase discovered can be thought of as a new 'universe', 
with its own rules of physics, and a corresponding set of 
'fundamental' excitations, or particles. At present there 
is great interest in electronic systems with ground states 
analogous to the nematic order found in liquid crystals. 
Order of this type is believed to exist in a wide variety 
of systems, including 2D frustrated antiferromagnetsi"— , 
the iron-pnictide materials^^— , Sr3Ru2 07i^, quantum Hall 
systems^ and there is strong evidence for quadrupolar or- 
der in '5f' actinide materials^. 

One particularly intriguing example is the quantum 
spin-nematic. In a spin-nematic, a spin-quadrupole mo- 
ment orders, breaking spin-rotation symmetry without 
breaking time-reversal symmetry^ ^^^d^" — . However, 
despite being first proposed more than 40 years ekgo^—, 
experimental evidence for the existence of this state re- 
mains scarce. One reason is that the common probes 
of magnetic order, including elastic neutron scattering, 
the static splitting of NMR lines and /xSR, measure the 
static internal magnetic field induced by magnetic or- 
der. These techniques are therefore blind to any state 
where time-reversal symmetry is preserved. This prob- 
lem is exacerbated by the fact that spin-nematic order is 
most easily stabilised at low temperature and high mag- 
netic field^, where many experiments are impossible to 
perform. Taken together, these facts represent a serious 
challenge - how can we probe for spin-nematic order in 
experiment? 

In this Letter we develop a general, quantitative, the- 
ory of NMR 1 /Ti relaxation rates in a quantum spin ne- 
matic, and show how it can be used to distinguish them 
from other, more conventional forms of magnetic order, 
even at high values of magnetic field. We find that, rel- 
ative to a conventional magnet, 1/Ti is (1) strongly sup- 
pressed, (2) does not exhibit a critical divergence ap- 
proaching a continuous phase transition [see Fig.[T], and 




FIG. 1: (Color online). a) 2-sublattice, bond-centered spin- 
nematic order with finite magnetisation, proposed to be the ground 
state of the quasi-lD magnet LiCuV04 for magnetic field h fa 43T 
[Tol4l2ll23l . The probability distribution for the moment on each 
bond is shown as a blue surface. The orthogonal 'directors' associ- 
ated with this type of nematic order are shown as red cylinders, b) 
Prediction of this Letter for the NMR 1/Ti relaxation rate in such 
a spin nematic in the region of its critical temperature T = Tpf , for 
parameters relevant to LiCuV04 [T^ . 



(3) exhibits unconventional power-law behaviour at low 
temperature. The analysis leading to these conclusions is 
based on symmetry considerations, and so remains valid 
regardless of the details of the system studied. To facil- 
itate comparison with experiment, all results arc stated 
in absolute units. 
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In recent years, there has been considerable progress 
in the understanding of spin-nematic states. A rich vein 
of candidate systems has been identified among spin- 1/2 
magnets in which ferromagnetic and antiferromagnetic 
interactions compete. Bond-centred spin-nematic order 
is found in the spin-1/2 Heisenberg model relevant to a 
family of square-lattice vanadatesi^r^^, and the multiple- 
spin-cxchangc model describing thin films of -^HeS^^. 
Also, the quasi-2D magnet NiGa2S4 has inspired the 
study of two different types of spin-nematic order in the 
spin-1 bilinear-biquadratic model on the triangular lat- 
ticei^"— . Perhaps the most promising candidate is the 
quasi-lD spin-chain material LiCuV04, which undergoes 
an unexpected low-temperature phase transition in mag- 
netic fields h « 43 ^10,11 ^ This has been explained in 
terms of the condensation of bound pairs of magnons 
to form 2-sublattice, bond-centred, spin-nematic order— 
[cf. Fig. Ufa)]. 

The high values of magnetic field needed to stabilise 
spin-nematic order in LiCuV04 render most experiments 
impractical. However measuring the NMR 1/Ti relax- 
ation rates of, e.g. "^Li, is technically feasibleii. The re- 
laxation rate can be expressed in terms of the magnet's 
dynamical spin susceptibility x"'^(q, w) js^ - ls^ ]. with the 
most general formic, 
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where 77V is the nuclear gyromagnetic ratio, wn is the 
NMR frequency, hext is the externally applied magnetic 
field, a and /3 label spin components, and J^"'^(q, hgxt) is 
a form factor describing the coupling between nuclear and 
electronic spins. It is therefore possible to look for the 
'finger print' of spin-nematic order through the coupling 
of its spin-dipole excitations to nuclear moments. 

In what follows we develop a quantitative, symmetry- 
based theory for the spin excitations of a quantum 
spin ncmatic in its critical region (T « Tjv), and low- 
temperature ordered phase (T ^ Iat), and use this to 
make corresponding predictions for 1/Ti. For concrete- 
ness, we consider a 2-sublattice 'antiferroquadrupolar' 
(AFQ) spin-nematic state in applied magnetic field, of 
the type proposed to occur in LiCuV04. We contrast 
this with the behaviour expected of a conventional canted 
antiferromagnetic (AFM). The results obtained are very 
general, and can easily be extended to quantum spin ne- 
matics in the absence of magnetic field. We note that a 
detailed theory for 1/Ti relaxation rates in LiCuV04, at 
higher temperatures where ID fluctuations predominate, 
has been developed in Ref. 36ll37} . 



We consider first both the overall scale of 1/Ti, and its 
behaviour in the vicinity of the critical point, T = T/y. At 
any continuous phase transition, the correlation length 
associated with the order parameter diverges as ^ oc t"'', 
where t = (T — Tn ) /Tn . In a conventional antifcrromag- 
net (AFM), the order parameter is a dipolc moment 
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FIG. 2: (Color online). Comparison of the NMR 1/Ti relaxation 
rate in a conventional antiferromagnet (AFM - red shaded area) and 
an antiferroquadrupolar (AFQ - blue shaded area) spin-nematic 
state. In the AFM, 1/Ti is seen to diverge approaching a critical 
point at T = Tm. Within mean field theory, 1/Ti cx J, where ^ is 
the correlation length. Predictions are taken from Eq. I IA6I I. with 
parameters relevant to KNiFs H: = 253X, T = i [T//iB]^, 
a = I, u = 1/3, K = 1/9, x±_ = 1/4 and F = 0.5. The re- 
sult has been convoluted with a Gaussian of standard deviation 
O.OSTm = 7.6X. In the AFQ, predictions arc taken from Eq. Illlt . 
and the parameters are the same as the AFM, with additionally 
A = 5 and 7 = 5. The 1/Ti relaxation rate is highly suppressed in 
the AFQ relative to the AFM and there is no critical divergence 
at T = Tm. The inset shows that there is instead a small bump 
at T = Tfj. The Feynman diagram describes the most important 
contribution to l/Ti relaxation rate [cf. Eq. l|10|l ] . 



of spin, and this divergence is accompanied by a criti- 
cal slowing in spin fluctuations. This in turn leads to 
a critical divergence in 1/Ti. Calculating this within 
time-dependent Ginzburg-Landau theorj*^"— , following 
Ref. we find (see supplementary materials). 
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where the form factor J^"'^(q, hgxt) — > -7^ is assumed 
to be constant, F is a phenomenological damp- 
ing parameter, K is a generalised elastic constant, 
f{x) ~ [arctanx — xl(\ + a:;^)]/2, A ^ tt is a momentum 
cut-off, lengthscales are measured in units of the lattice 
spacing, energy in units of /cbTn and time in units of 
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h/{kQTu). Approaching the critical point, Eq. (jA6[) pre- 
dicts that 1/Ti diverges as 1/Ti cx ^ (xt^" . 

More generally, the relaxation rate diverges as^. 



z+2-d ^ :^{d-2-z)v 



(3) 



where, within mean field theory, the dynamical exponent 
z = 2 and correlation length exponent v = 1/2. Ex- 
periment confirms that, at a critical point separating an 
AFM fi"om the high-temperature paramagnet, 1/Ti is in- 
deed described by Eq. ^ 

In order to understand the role of critical fiuctuations 
in a spin nematic, it is necessary to generalise this anal- 
ysis to a quadrupolar order parameter. A continuous 
phase transition from the paramagnet into the AFQ state 
is allowed by symmetry, and we find that the energy cost 
of static fluctuations of AFQ order is given byi^i^. 
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where F is a canting field describing spin-dipole fluctu- 
ations parallel to applied magnetic field, an irrelevant 
coupling hF has been ignored, and the order parameter 
for the 2-sublattice AFQ state is given by. 



Q±(r) = 



Q^/(r)-Q^3^(r) 



(5) 



where = (Sf)^ - {Sl)\ = SfSj + 5^5? and 

z — A, B. In order to separate longitudinal and trans- 
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verse fluctuations we write Ql ~^ ~ ?y[(ri,J)^ — (n^)^] and 
= 2irjn^ri[, where = (nj,n^,0) is a unit vector 
aligned with the nematic directors and ha • ne = 0. This 
results in, 
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with the director stiffness, Pii{rf) = rfK. 

The longitudinal fiuctuations described by r]{r,t) do 
not couple directly to the nuclear spin lattice because 
the order parameter Q_L(r) is time-reversal symmetric. 
The 1/Ti relaxation rate is instead controlled by trans- 
verse fiuctuations of the order parameter, parameterised 
by rii and P(r, t), which arc accompanied by a small, fiuc- 
tuating, dipole momenti^. Writing Ha = {^/l--(j)^,(f>,0) 
and hb = (—0, \/l — (f>^, 0) and linearising the resulting 
equations of motion, we find. 
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9t/^(r,t)« ^V>(r,t) + ^V^P(r,t) + 0(r,t), (7) 
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where F, 7 and A set the relaxation timescales for the 
fields ?7, (j) and P, and the mean field approximation 
w —> has been assumed. The auxiliary fields C??, C(j> ^-nd 
C;, which parameterise the damping of long- wavelength 
fiuctuations by interaction with short-wavelength excita- 
tions, obey the correlation functions, 

(C^(r, t)Ur',t')) = 2rfcBT<5(r - r')S{t - t') 
(C0(r, t)C^{r', t')) = 2-fksTSir - r')<5(t - t') 
Uv,tMr',t')) =2XkBTV'Sir-v')Sit-t'). (8) 

Fluctuations of the total magnetisation commute with 
'Hafq [Eq. dH)], so the q = component of P is conserved. 
Solving the equations of motion, Eq. (jBSp . results in a 
dispersion. 
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where the first term in describes the coherent motion 
of the quadrupolcs and the second term is dissipative. 

We are now in a position to calculate the dynamical 
spin susceptibility of the spin-nematic. This is not given 
by the dynamical correlations of the order parameter, 
but rather the associated fiuctuations of the spin-dipole 
moments. 

6S^{r,t) oc Z^(r,i) + [(77)+,5?7(r,t)]P(r,t) w 6ri{v,t)l\v,t), 

where the (5S'^(r, t) cx P(r, t) terms are negligible for l/Ti 
relaxation. This is because 1/Ti probes the spin suscep- 
tibility at w — > [cf. Eq. H])], and at low energies the 
P term is suppressed by the fact that /^(q = 0,i) is a 
conserved quantity. The dynamical spin susceptibility is 
calculated from the diagram shown in Fig. [21 and this 
leads to, 



d^k f duj' 
{keTN)^ J J 27 



(10) 



where Xq = K~'^{^^'^ + q^)~^. Physically, these 
equations describe a longitudinal fluctuation of the 
quadrupole order parameter, which then cants, result- 
ing in a spin-dipole fluctuation. In order for this process 
to occur, it is necessary that 3> tj, where the correla- 
tion time for the P field is given by Ti(q) w 2x^/Aq^ in 
the paramagnet and T^(q) ~ r~^Xq. Close to the critical 
point this is equivalent to requiring x^FA'/A ^ 1. 

Substituting Eq. PH)) into Eq. ([T]) results in a relax- 
ation rate. 
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where the form factor, J^, is again assumed constant. 

A numerical integration of this expression for an 
isotropic, 3D AFQ state results in the temperature de- 
pendence shown in Fig. [21 In direct contrast to the AFM, 
there is no divergence at the critical point, but instead 
a small bump. The relaxation rate in the AFQ is sup- 
pressed relative to the AFM even away from the critical 
point. For t w 0{1), assuming comparable parameters, 
the ratio of the two is small as, 



1 



rpAFQ' yAFM 



< 1. 



(12) 



The enhancement of the relaxation rate in the paramag- 
net relative to the ordered phase (cf. Fig. [3]) is due to 
the greater phase space for fluctuations. 

We now turn to the low-temperature ordered state. 
In a conventional AFM, at low temperatures, l/Ti is 
driven by the Raman scattering of thermally excited 
spin waves^2iH. In the absence of magnetic anisotropy, 
a canted AFM in field has a single linearly-dispersing 
Goldstone mode, associated with rotations of the order 
parameter about the direction of magnetic field. This 
controls 1/Ti for T <C TV and leads, in 3D, to 
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where rus is the staggered ordered moment, V^eW is the 
volume of a unit cell, Vs = (^'a^'b^'c)^ is the mean spin- 
wave velocity and XJ- th^ transverse susceptibility. A 
constant form factor J- has been assumed. 

The NMR l/Ti relaxation rate at low temperatures 
in a quantum spin nematic is controlled by the long- 
wavelength excitations of the nematic order parameter. 
We have constructed a sigma-model model description of 
the low-energy excitations of AFQ order, based on the 
SU(3) algebra needed to describe general rotations of the 
quadrupolar order parameter. Details of this work will be 
presented elsewhere^ - here we concentrate on the impli- 
cations for NMR 1/Ti relaxation rates in a 2-sublattice 
AFQ state of the type proposed for LiCuV04 [lollTl| . 

In the 2-sublattice AFQ there are four long- wavelength 
modes, one of which is a linearly dispersing Goldstone 
mode associated with rotations of the quadrupole mo- 
ments within the plane perpendicular to magnetic field. 
There are also three gapped modes, two associated with 
out-of-plane rotations of quadrupole moments, and one 
that has the character of a fluctuating spin-density wave. 
This last mode has a strong dipole moment, but is gapped 
even in the absence of magnetic field. 



At temperatures small compared with all energy gaps, 
only the Goldstone mode contributes to relaxation. Spin- 
dipole fiuctuations are associated with time derivatives of 
the field, and, in 3D and for T <C T^r we find. 
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where Iht 

the susceptibility associated with the Goldstone mode, 
and the algebraic mean of its velocity^. Comparing 
Eq. and Eq. ([Ti]) . we see that the low-temperature 
prediction oc for a spin-nematic state in 

magnetic field is dramatically suppressed relative to the 



1/T( 



AFM 



predicted for a canted antiferromagnet. 



We note that relaxation due to excitation of acoustic 
phonons also has a oc T'' behaviour—, but could 

in principle be distinguished through its (lack of) mag- 
netic field dependence. 

The temperature dependence of l/Ti in the spin ne- 
matic crosses over to 1/Ti oc when the temperature 
becomes comparable with the energy gap, Az, to dynam- 
ical spin-density wave fluctuations, which couple more 
strongly to nuclear spinai^. We also also find 1 /Ti oc 
for a 3D AFQ state in the absence of magnetic field, in 
agreement with the 1/Ti oc T^'^'^ quoted in Ref. 
However the 1/Ti oc found in magnetic field arises 
from a different physical process (see supplementary ma- 
terial). 

In conclusion, we have shown that the NMR 1/Ti re- 
laxation rate can be used to distinguish a spin-nematic 
state from conventionally ordered magnets. We find that: 
(1) 1/Ti is strongly suppressed relative to conventional 
magnets; (2) unlike conventional magnets, l/Ti does not 
diverge at a continuous phase transition; and (3) l/Ti ex- 
hibits unusual power-law behaviour at low temperature. 
The intrinsic smallness of 1/Ti in quantum spin nemat- 
ics may render low-temperature power laws difficult to 
observe. However, the absence of a critical divergence 
in 1/Ti is a robust, qualitative feature, which should be 
experimentally testable even at the high magnetic fields 
and low temperatures needed to stabilise spin-nematic 
order in systems like LiCuV04. 
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Appendix A: Theory of NMR relaxation rates in a 
conventional canted antiferromagnet 

For completeness and ease of comparison, we develop 
a theory of the relaxation rate in a conventional canted 
antiferromagnet (AFM), where 1/Ti is known to diverge 
for T ^ Tn [ll-til- Spin excitations in the AFM can 
be described phenomenologically using time-dependent 
Ginzburg-Landau theory^^— , written in terms of the order 
parameter ms(r), and its associated canting field P(r), 
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at 9 
2 ' 
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(Vms 
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z\2 



2X± 



{I 
(Al) 



where a term hl^ is ignored as irrelevant. Here 
i = {T — T[^)/T!^, K is a. generalised elastic constant, x± 
the transverse susceptibility, a > and u > are posi- 
tive constants and magnetic field is applied parallel to the 
z-axis. Length scales are measured in units of the lattice 
spacing and energy scales in units of ksT^. Expanding, 

ms(r) = 77(r)n(r) = [(77) + S'n{r)]n{r), 

where n(r) is a 2D unit vector, results in, 
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where the spin stiffness is defined by PsiTf) ~ Krf'. The 
Ginzburg-Landau theory predicts a diverging correlation 
length ^ cx t~'^ approaching the critical point, with mean 
field exponent v = 1/2. 

The critical divergence of ^ is accompanied by a criti- 
cal slowing down in fiuctuations of AFM order. This can 
be described by adding dynamics to the static Ginzburg- 
Landau model, to capture the way in which the slow, 
long- wavelength fluctuations of AFM order are damped 
by interaction with short-wavelength fluctuations. Fol- 
lowing Ref.'s [B-ll], we write, 

dtS-niv, t) = -TK{r^ - V')S7j{r, t) + C(r, t), (A3) 

where time is measured in units of h/keTj^, F^^ sets the 
rate of damping, m — > 0, and C(i"ii) is a- (white) noise 
term with correlation function. 



(C(r, t)C(r', t')) = 2fcBrr,5(r - r')S{t - t'). 



(A4) 



It follows that the correlation time is given by 
T,,(q) « F-^Xq, with xq = K-^{£,-^ + q^)"^ Approach- 
ing Tn, t^(0) diverges as t^(0) cx with dynamical ex- 
ponent z = 2. 

In the case of the antiferromagnet, fluctuations of the 
field T] are directly related to fiuctuations of the spin- 
dipole moments, since, for moments canted in the x- 
z plane, 5^(r,t) cx ±77(r, i). Therefore, the dynamic 



susceptibility associated with r) controls the critical be- 
haviour of the 1/Ti relaxation rate. The fluctuation- 
dissipation theorem allows the imaginary part of the dy- 
namic susceptibility to be calculated from Eq. (jA3p and 
Eq. (gH) as. 



hu! 



(A5) 



where g\ is the Landc g-factor. Thus the dynamic suscep- 
tibility diverges approaching the critical point for both 
q ^ and a; — ?> 0. 

The mean- field behaviour of the 1 /Ti relaxation rate 
can be calculated from Eq. ([l} in the main text and 
Eq. (|X5l) as. 
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(giMB)^ _,2x- 



ksT 
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(ksT^r FA'2 



^f{Am, (A6) 



where a constant form factor, J^"'^(q, hgxt) ^ -7^, 
is assumed, A is a momentum cut-off and 
/(x) = [arctana; — a;/(l -I- a;^)]/2, where f{oo) = 7r/4. 
Approaching the critical point Eq. (|A6|) predicts 
1/Ti cx ^ cx i-", i.e. the NMR l/Ti relaxation rate 
diverges with mean field exponent ly ~ 1/2. More 
generally this can be written as^, 



— cx oc t 



(d-2-z)v 



(A7) 



where the exponents z and v exponents take on values 
appropriate for an 0(2) phase transition. However ex- 
periment confirms that mean field theory is qualitatively 
correct in predicting a divergence in l/Ti [2l-j3- 



Appendix B: NMR 1/Ti in LiCuV04 at high 
magnetic field 

It has been suggested that a 2-sublattice antiferro- 
quadrupolar (AFQ) spin nematic state may be realised 
in LiCuV04 at high magnetic fieldiSiii [see Fig. 1 in the 
main text]. We here extend the analysis of 1/Ti in the 
AFQ state presented in the main text to the specific ex- 
ample of LiCuV04. The difference in the analysis is that 
LiCuV04 is quasi- 1-dimensional, and, in consequence, 
the stiffness, correlation length and all related proper- 
ties are anisotropic. 

The anisotropic Ginzburg-Landau Hamiltonian [cf. 
Eq. ([¥]) in the main text] is given by. 



H 



-Qi 
2 



K, 



-Ql 



1 



in' 



(Bl) 



where ji = {a,b,c} labels the crystallographic 
axes, the Einstein summation convention is assumed, 
i = {T — T|\i)/rN is the reduced temperature, u > 0, Xz 



7 



is the perpendicular susceptibility and z defines the field 
direction. 

The order parameter can be rewritten using 

<3f = ^[("?)' - and Qf^ = 2777ijn^, where 

rij — (n^, n1, 0) is a unit vector aligned with the nematic 
directors and ha • ne = 0. This results in, 



n = 



at 2 , Kf^ 



(B2) 

where the director stiffness, p^{ri) = rfK^, is spatially 
anisotropic. 

The phenomenological equation of motion associated 
with the longitudinal order parameter field is given by, 

dtSiiiv, t) = -VK^{i-^ - dl)5ri{v, t) + C(r, t), (B3) 



where 



K^/ai in the paramagnet, 
~ \l ^K^j2at in the ordered state and, 

(C(r, mv',t')) = 2ksTT6ir - r')(5(t - t'). (B4) 
The equations of motion for the transverse fields are, 
P(r,t) 



X: 



Q 



dtnv.t) « p^a>(r,t) + -|9^P(r,i) + 0(r,i), (B5) 

where is anisotropic and, 

(C^r, t)C4r', t')) ^ 2jkBTS{r - r')5{t - t') 
(Ci(r,t)Ci(r',i')) = 2X^kBTdlSir - r')<5(t - t'). (B6) 



Solving these equations results in the dispersion relation, 

(B7) 



2„2 ^^/^ 



with, 



A:, co' 




k+q, co+co' 



FIG. 3: (Color online). Feynman diagram used to calculate the 
dynamic spin susceptibility close to the critical point in the anti- 
ferroquadrupolar (AFQ) spin nematic state. A red line is used to 
denote the order parameter field Srj, while a black line denotes the 
canting field P. 



The NMR relaxation rate can be calculated from the 
Feynman diagram shown in Fig. [31 and is given by. 
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2 {gm? {ksTy 
7n o 



(/cbTn 



dkidk2duj 
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where Xki K ^{^ ^ + hf) ^ and the approximation 
of a constant form factor, J^"''(q, hgxt) — > J^, is justified 
for '^Li NMR in LiCuV04 in Secti on Q b elow. Length- 
scales have been rescaled using y^Kf^/K, and therefore 
the momentum integral can be taken over a sphere. The 
expression for l/Ti given in Eq. (|B10[) is equivalent to 
Eq. ((TT|) in the main text, except for the fact that the 
correlation length, ^, stiffness p^, damping A and all re- 
lated quantities have been replaced with their geometric 
mean. 

The result of numerically integrating this equation for 
parameters appropriate to LiCuV04 is shown in Fig. |4l 
where the response has been convoluted with a Gaussian 
of width O.OSTn. This clearly shows that there is no crit- 
ical divergence at T = Tn, merely a small bump. For a 
reasonable set of parameters (quoted in the figure cap- 
tion) the 1/Ti is also extremely small compared with a 
conventional AFM. 



and, 



(B8) 



Appendix C: The form factor for ^Li NMR in the 
spin-nematic state 
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where K — {K^Ki^K^Y/'^ and p ~ Krf are geomet- 
ric means and we have made the assumption that 
A,, = \K^/K. 



In Section [B1 above, we assumed the NMR form factor 
to be a constant. We now justify this assumption for the 
case of ^Li NMR in LiCuV04. 

We assume that the Li nuclear spin interacts via a 
transferred hyperfine interaction with the electronic spins 
on the 4 nearest- neighbour Cu atoms [see Fig. [5] . The 
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FIG. 4: (Color online). Predictions for the NMR l/Ti relaxation 
rate in the spin nematic phase proposed to exist in LiCuV04 at 
high magnetic fieldifliii. There is no critical divergence at T = Tm , 
but instead a small bump. Predictions are taken from Eq. IIBlOl l 
with Tm =2K,r = 2.5 x IQ-^ [T//^b]^, a = 1, u = 1/3, K = 0.22, 
XX = 1/20, r = 0.5, A = 5 and 7 = 5, where the static parameters 
have been adapted from Ref . Il2f . The result has been convoluted 
with a Gaussian of standard deviation O.OSTm. 



internal magnetic field at the Li site is given by, 

i 

where is the nuclear-electron coupling tensor and 
Ta.i{t) is the magnetic moment associated with the zth 
Cu atom in the 4-site plaquettc. The form factor for 
1/Ti relaxation is given byi^, 



7,5 



AlOL A 



5/3 

-q' 
(C2) 



where Rhe^t is a rotation matrix relating the direction of 
the external magnetic field to the crystallographic coor- 
dinate axes and, 



al3 



(C3) 



In the spin nematic phase the Cu spin-l/2's form 
time-reversal-symmetric triplets, and the quadrupolar or- 
der parameter is associated with the bonds between Cu 
atoms. Excitations of this state mix a small spin-dipole 
component into the wavefunction, and these spin-dipole 
fluctuations are also associated with the bonds. The in- 
ternal magnetic field at the Li site can be expressed in 
terms of the bond spin fluctuations by, 



(C4) 



where A13 = [Ai + ^3)72, ^24 = {A2 + Ai^jl and 
ijj(t) describes the fluctuating spin-dipole moment as- 



m, 



sociated with the bond joining the site i with the site j 
[see Fig. [S]for the numbering of Cu atoms]. Using the 




FIG. 5: (Color online). The local environment of a Li atom in 
LiCuV04. The internal magnetic field, /ijnt, at the Li site is as- 
sumed to arise predominantly from a transferred hypefine interac- 
tion with the four nearest-neighbour Cu atoms. In the antiferro- 
quadrupolar (AFQ) spin nematic state, the quadrupolar order pa- 
rameter exists on the bonds between neighbouring Cu atoms [see 
Fig. 1 in the main text] . Fluctuation of the order parameter creates 
spin-dipole fluctuations on these bonds. These drive fluctuations 
of /lint the Li site, and thus couple to nuclear 1/Ti relaxation. 



symmetry of the Li environment, these tensors can be 
expressed as. 



t-13 



A, 



•24 
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and in consequence, 



where. 




2 



(C6) 



(C7) 



The dominant spin-dipole fluctuations are in-phase on 
neighbouring bonds (ie. q « 0) and parallel to the di- 
rection of applied magnetic field. In consequence, the 
relevant form factors are given by, 

\2 







(C8) 



Thus, for magnetic field applied in the a and c 
directions the form factor is well approximated by 
a constant. The nuclear-electron tensor can be 
probed in the low-field, magnetically-ordered phase 
to find A?"^ w 0.025r/^B consequence 
J""(q w 0, h^) w J"==(q w 0, h^) w 2.5 X 10-3[T//ib]^. 



Appendix D: NMR 1/Ti at low temperature in the 
spin nematic state 

Low temperature measurement of the NMR 1/Ti re- 
laxation rate could also be used to distinguish the spin- 



9 



nematic state from conventional magnets. At low tem- 
peratures 1 /Ti is activated above an energy gap and then 
follows a power law behaviour as a function of T , with 
different exponents in the spin nematic compared to con- 
ventional magnets. We consider here the spin nematic 
state that may be realised in LiCuV04. 

In a conventional coUinear antiferromagnet at low 
temperaturi 



,13 



1 



kp.T 



(Dl) 



where a constant form factor J- is assumed, is the 
ordered moment, V^eii = aofcoco is the volume of a unit 
cell, A is an energy gap, Vs is the geometric mean of the 
spin wave velocities, X-L is the perpendicular susceptibil- 
ity and. 



$(a;) = .T2Lii(e-i/^) + x^LizCe-^/^^), 



(D2) 



with Li„i(z) — X^i^o^V'™ tiic m*'' polylogarithm of z. 
For temperatures T S> A, 



Ti ^ 487r 



2 2 iW 



while for temperatures T ^ A, 



cell 



(D3) 



(D4) 



For a discussion of cases where the form factor is not 
constant see Ref. [Tsj . 

In a 2-sublattice AFQ state in large applied magnetic 
field, 1/Ti can be calculated by an analogous method. 
The excitations of the quadrupolar order parameter can 
be described by a non-linear sigma model field theoryi^, 
and this can be linearised to give the effective action. 



1 



/3 



■^Vcell Jo 



of the order parameter in the plane perpendicular to the 
magnetic field, and these are described by the field ipi. 
The gap Agx allows for a small exchange aniso trop y in 
this plane, as has been measured in LiCuV04 ^M,- At 
higher temperatures three additional modes with much 
larger energy gaps become important. The fields (j}Q and 
(/itt describe in-phase and anti-phase rotations of the or- 
der parameter towards the direction of the magnetic field. 
The field 04 describes fluctuations of a spin-dipole mo- 
ment parallel to the field direction. It is not clear with- 
out solving a specific microscopic model, which of the 
gaps, Az(/i), Axy,o(^) or Axy,7r(/i), is smallest. However, 
as far as the 1/Ti relaxation rate is concerned, the three 
modes behave in an equivalent manner for ^Li NMR in 
LiCuV04. 






k^+k2+q 



FIG. 6: Diagrams used to calculate the low-temperature, dynamic 
spin susceptibility in the spin-nematic phase proposed to exist in 
LiCuV04. At low temperatures the dynamic spin susceptibility is 
controlled by the field and is calculated using diagrams (a), (b) 
and (c). Slashes denote differentiation of the fields with respect to 
time. It follows that 1/Ti oc T^. At temperatures comparable to 
the gap Az, excitations described by the field <j)4 become important, 
and diagram (d) is used to calculate the dynamic spin susceptibility. 
It follows from this diagram that 1/Ti oc T^. 



2^1,2 
4 



Xxy,7r(/i)(9r07r)^ 



P.y,.Ah){df.<l>.f + x.yAh)A%A'')€}, (D5) 



The action can be used to calculate the dynamic spin 
susceptibility, and l/Ti follows from Eq. ([1]) in the main 
text. At low temperature we just consider the field 0i. 
The lowest order contribution to the relaxation rate is 
the 'three magnon' term, which involves averages over 
six (jji fields with two time derivatives. Evaluating this 
gives. 



where p's denote order parameter stiffness', x's suscep- 
tibilities, A's energy gaps and /.j = {a,6,c}. This action 
describes four gapped modes with dispersia of the form 

Wk = \J + v^kf^, where = ^IhJx- At low tem- 
peratures physical properties are controlled by rotations 
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where (mf)^ = V3{g^lB? m±}\ /2 and, 



where, 



A! 



J dxidx2]lx-f 
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(D7) 

where = (v?a'^?b''^?c)^^'^ is the geometric mean of 



the velocites and = y p^^/x?- At temperatures 
T <^ Aex this can be approximated by, 

Ti 1287r5 V s ^ w ^Q(^Q)9 ' ^ ^ 



where. 



^(^)-w(^) Jl^dx.dxJxi 



2 / 



fcBT 



k^T 



{xi + 0:2)2 - 
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■ . X, X, , r- (D13) 

(e^i - l)(e^2 _ i)((.xt+x^ V / 

At temperatures T ^ A^ this can be approximated by, 



where, 



Jo 



For temperatures T ^ Aex, 
J- 



1 ^ r^effx2^2 (fceT) 
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where, 



For temperatures T Az 



/hT = / dxidX2 XlX2{xi + X2) [Xi + X1X2 + X2) 

Jo 

pXi+X2 

X — — ■ TV ~ 146.5. 



(e^i - l)(e^2 _ l)(e^i+^2 _ i) 



(Dll) 



At higher temperatures the gapped modes associated 
with the fields 04, (j)Q and c/)^ become important, and 
dominate the 1/Ti relaxation. It is not a priori clear 
what the relative size of the three gaps is, or how they 
compare to Tm. However, the contribution to the 1/Ti 
relaxation rate associated with each mode is the same, 
and so we just consider in detail the ^4 mode. The con- 
tribution to the relaxation rate is given by. 



Ti IGttS 
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where. 



IhT = / dxidx2 XiX2{xi + X2) 

Ja 



.66. 

(D17) 



We have shown that the NMR l/Ti relaxation rate 
has a distinctive temperature dependence at low T in the 
AFQ spin ncmatic state. In principle this could be used 
to recognise the spin ncmatic state. However, this may 
prove experimentally difficult due to the long relaxation 
times. 
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